Spin-texture and magnetic anisotropy of Co adsorbed Bi2Se3 topological insulator 

surfaces 
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Based upon first-principles methods, we investigate the magnetic anisotropy and the spin-texture 
of Co adatoms embedded in the topmost Se network of the topological insulator Bi2Se3 surface. We 
find the formation of energetically stable magnetic moment perpendicular to the surface plane, Sz. 
Our results for the pristine Bi2Se3 surface indicate the presence of helical spin-texture not only in 
the massless surface Dirac states, but also surface states resonant within the valence band present 
spin-texture. On the other hand, upon the presence of Co adatoms we find that the out-of-plane 
surface magnetism represents the dominant spin state (Sz), while the planar spin components, Sx 
and Sy, are almost suppressed. 

PACS numbers: 73.20.At, 75.30.Hx, 85.75.-d, 71.20.-b 

Topological insulators (TIs) compose a new quantum 
matter phase in solid state physics that may provide 
an emplacement for fundamental physics understanding 
as well bases for novel technologies, such as topological 
quantum computing, information processing and spin- 
tronic applications. Recently 3D TIs were predicted by 
Fu et al. [l| and experimentally observed by Hsieh et al. 
[3 in Bii-j^Sbj:. Further, other compounds like Bi2Te3, 
BiaSes, SbzTeg, and TlBiSca were identified as 3D TIs. 
[SHOJ In these materials the spin-orbit (SO) interaction 
rules the formation of topological states. While the bulk 
is insulating, on the edges there are robust conducting 
states. Indeed, on the surface they present a (single) 
Dirac crossing at the F point. They show a spin-texture 
that allows spin currents without dissipation, supress- 
ing back-scattering if time reversal symmetry (TRS) is 
preserved. In this way the understanding of how the 
protected surface states behaves in the presence of impu- 
rities is an important subject to be explored, especially 
for technological applications. Magnetic impurities can 
destroy the robust surface metalicity by breaking TRS, 
while nonmagnetic impurities maintain TRS. The inter- 
action between a magnetic moment and the topologi- 
cally nontrivial surface states opens a surface band gap. 
[lot lllj In this context an important issue to be inves- 
tigated is the energetic stability of the spin orientation, 
and the changes on the spin-texture upon the adsorption 
of transition metal on TI surfaces, which is especially 
desired for spintronic and information processing appli- 
cations. 

In this letter, through the calculation of magnetic 
anisotropy energy (MAE) we find that Co adatoms ly- 
ing on the TI Bi2Se3 surface exhibit energetically sta- 
ble magnetic moment perpendicular to the surface. Fur- 
ther, spin-texture calculations reveal that the Co adatom 
exhibits spin component ^100% polarized out-of-plane. 



whereas the helical spin-texture of the massive Dirac 
cone is drastically reduced. In addition, in consonance 
with previous experimental findings, the presence of Co 
adatoms breaks TRS, opening up a surface band gap, 
giving rise to massive Dirac fermions. 

Our results were obtained performing ab-initio elec- 
tronic structure calculations, using the density func- 
tional theory (DFT) within the generalized gradient 
approximation (GGA) for the exchange and correla- 
tion potential. [12] The SO interaction have been self- 
consistently treated by using fully relativistic pseudopo- 
tential within the projector augmented wave method 
(PAW).[13| We use the Vienna Ab initio Package Simu- 
lation (VASP).il4. A plane wave basis set is used, with 
a cut-off energy of 300 eV. The Brillouin-zone is sam- 
pled, according to the Monkhorst-Pack[15J scheme, by 
using a number of k-points such that the total energy is 
converged. The (111) surface of Bi2Se3 was investigated 
using slab method with a vacuum layer of at least 8A. 
Here we have used a thickness of 4 quintuple layers, keep- 
ing the experimental lattice parameter. We verify that 
this slab is enough to obtain the correct surface massless 
Dirac cone. 

In Fig. [T]we show the band structure around the bulk 
band gap, and the surface spin-texture when a Co atom is 
adsorbed onto Bi2Se3 surface. The surface spin-texture is 
affected by the presence of the magnetic impurity, as we 
can see in Fig. [Ij-b-j. The original massless Dirac cone 
(black dotted lines) is now splitted (purple line) where 
the lower massive Dirac cone lies below the bulk valence 
band, and the upper massive Dirac cone stays inside the 
bulk band gap. Fig. [T^. Thus, indicating that the TRS 
has been removed upon the presence of Co adatoms on 
the TI surface. 

The calculated expectation values of Sx and Sy of the 
surface protected massless Dirac cone for the pristine 



Bi2Se3 (black lines of Fig. [T|) are ~60% spin-polarized 
with respect to ?i/2, in agreement with a recent calcula- 
tion. [16] By looking at Figs. [T]-e and 1-f we see that the 
Sx and Sy for the electronic states of the upper massive 
Dirac cone far from the F point of the Brillouin zone, are 
somewhat the same as compared to the ones of pristine 
Bi2Se3 surface. On the other hand, for wave vectors near 
the F point, fc — > F, those spin components tend to zero. 
That is, the projection of the spin-components on the 
xy plane becomes negligible. Meanwhile, as depicted in 
Fig-HJ-g, there is an increasing on the out of the plane, Sz, 
component close to the center of the Brillouin zone, which 
can be attributed to the formation of Co-Se chemical 
bonds. Here the expectation values of spin components 
are obtained from {Sa{k)) = h/2{^{k)\aa\^{k)), where 
"^(k) is the spinor wavefunction and CTq are the Pauli ma- 
trices. These results are in agreement with recent theo- 
retical calculations based on few band phenomenological 
models. ^^ 

We find (five) d orbitals inside the bulk band gap 
mostly localized in the Co atom (Fig. [T]-a) . We can sepa- 
rate these orbitals in two classes according to their sym- 
metry. Two of them (red and pink lines) suffer a strong 
hybridization with the surface Se atoms, resulting in a 
quasi-Dirac crossing at the F point. These two orbitals 
have Se-pz and Co-dxy and -d^^-y^ symmetry characters. 
As we can see from Fig. [T]-h and [T]-i they present small 
Sx and Sy spin components, and the left-handed spin he- 
licity has been maintained, while there is an out of plane 
Sz of around h/A for the electronic states near the edge 
of the Brillouin zone. Fig. [T]-j. The other three Co-Sd 
orbitals (green lines of the Fig. [T]-a) are more localized, 
being two occupied bands {dxy-yz) and one unoccupied 
{dz^). The magnitude of the Sz component of these three 
bands are ^ h/2 spin aligned perpendicular to the TI sur- 
face (see Fig.[T]-d), and do not dependent on k-direction. 
In contrast the Co-dxy and -dx^-y^ (Fig. [T]-j), due to the 
hybridization, the spin polarization is dependent on the 
k-direction. 

It is worth to point out that the results presented in 
Fig. [T] are for a Co coverage of 1/9 monolayer (ML). 
In addition we have considered coverages of 1, 1/4 and 
1/16 ML, as schematicly depicted in Fig. EI For 1/9 ML 
Co coverage on the surface of Bi2Se3 there is not new 
Dirac point, but two bands (red and pink lines of this 
figure) are getting close at the F point. By increasing 
the Co coverage to 1 or 1/4 ML, this two bands are split- 
ted. However by decreasing the Co coverage to 1/16 ML 
the two bands touch to each other at the F point, forming 
a new Dirac cone, similarly as experimentally observed 
for Fe deposited on Bi2Se3. (201 We verify that the spin- 
texture for the other bands (then the red and pink ones 
from Fig. [ij are similar for all Co coverages. 

In order to understand the spin projection as a function 
of the energy states we plot in Fig. [3] the Spin-Density 
projected for each spin component (a — x,y,z). This 
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FIG. 1: (Color online) (a) Band structure around the bulk 
band gap when a Co atom is adsorbed on the surface of the 
Bi2Se3. The upper and lower shadow regions are the bulk 
conduction band and bulk valence band. The black dotted 
lines are the original massless Dirac cones without the pres- 
ence of the Co adsorbed, and the purple line is the splitted 
massive Dirac cone. The red, pink and the three green lines 
are the bands introduced by the Co impurity, (b to J) Ex- 
pectation values of the spin components {Sa{k)) (a — x, y, z) 
along the K^j-F-My direction, (b-d) the spin of the three up- 
per Co-d orbitals green bands of (a); (e-g) the massive Dirac 
cone; and (h-j) are the lower Co-d orbitals (pink and red from 
(a)). Also is plotted (thin black lines) the expectation values 
for the pristine Bi2Se3, where the dotted black lines are the 
components for the (unoccupied) upper massless Dirac cone 
and the full black lines are the components for the (occupied) 
lower massless Dirac cone. At the right-bottom corner are 
the calculated K-F and F-M directions represented in the 
first Brillouin zone. 



projection was computed for a mesh of k-points along 
the Kx-T-My direction. For pristine Bi2Se3 (Fig. [3]- 
b and -c) we observe that the occupied massless Dirac 
cone presents higher density of spins than the unoccu- 
pied one. From this figure we also observe that there is a 
spin-texture for states inside the bulk valence band reach- 
ing '-^l eV below the Fermi level. These spin-textures 
are due to the presence of surface states resonant within 
the bulk valence band. When the Co atom is adsorbed 
on the TI surface (Fig. |3]-d and -e) the {Sx) and (5*^) 
components reduce drastically, while the {Sz) increases, 
showing that the out-of-plane surface magnetism are the 
dominant spin states. 

The Co atom adsorbed on the Bi2Se3 surface is bonded 
to three top surface Se and three Bi atoms of the second 
layer, as showed inset of Fig. 01 This system is exother- 
mic by 2.09 eV per Co atom with respect to clean surface 
and isolated Co dimers. There is an attractive inter ac- 
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FIG. 2: (Color online) Schematic representantion of Co ad- 
sorbed on the TI surface for 1, 1/4, 1/9 and 1/16 monolayer 
coverages, from (a) to (d), respectively. The whiter balls rep- 
resent the Se topmost surface layer, and the red balls represent 
the embedded Co atoms. 
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FIG. 3: (Color online) Band structure (a) and Spin-Density 
projected on each So- spin component for pristine Bi2Se3 along 
the r-K (b) and F-M (c) directions, (d-e) The Spin-Density 
projections for the system when a Co atom is adsorbed on 
the TI surface. The positive and negative signs for <So- > are 
with respect to the K^^-F-Mj, directions, as indicated inset at 
the hexagon diagram. The left-handed helicity shown in the 
hexagon diagram is for a cut above the Fermi level. 
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FIG. 4: (Color online) Adsorption energy as a function of 
Co- Co average distance distributed on the surface of Bi2Se3. 
The adsorption energy is computed with respect to the clean 
surface and isolated C02 dimer. Inset is the atomic arrange- 
ment when a Co atom is adsorbed on the Bi2Se3 surface. The 
Co atom is bonded to the first (Se) and second (Bi) atomic 
layers. 



can works as a local magnetic field inducing an anoma- 



lous Hall effect. [HE 



The spin-texture of the TI~magnetic-impurity system 
is dependent on the spin polarization. In order to ver- 
ify the energetic stability of our calculated spin-texture 
(Fig. [T] and |3]) we calculate the magnetic anisotropy en- 
ergy of the magnetization induced by the Co adatoms 
on the Bi2Se3 surface. As we can see in Fig. [5] the min- 
imum energy occurs for the magnetic moment aligned 
perpendicular to the TI surface, and the energy barrier 
is verified when the magnetic moment is in-plane. A sin- 
gle Co impurity adsorbed on TI surface presents a MAE 
of 6 meV per impurity. 
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tion between the Co atoms when they are adsorbed on 
the Bi2Se3 surface. By increasing the Co coverage, as we 
can see in Fig. 2] the adsorption energy increases. The 
Co-Co interaction favors a FM state. The total energy 
difference per pair of Co atoms with spin parallel as com- 
pared to that of spin anti-parallel is 57 meV favoring a 
FM configuration. A FM phase was theoretically pre- 
dicted for magnetic inipurities in TI when Sz spin com- 
ponent is dominant. [ITl ll9| This out-of-plane spin net 



FIG. 5: Magnetic anisotropy energy computed as the total 
energy as a function of the magnetic moment alignment. The 
zero angle is for a magnetic moment aligned out-of-plane, per- 
pendicular to the TI surface (z direction), xy is the basal TI 
plane. 

In summary our first principles calculations provide a 
fundamental understanding of TI properties upon the ad- 
sorption of a magnetic impurity. The adsorbed transition 
metal Co impurity, give rise to massive Dirac fermions. 



inducing an energetically stable out-of-plane spin com- 
ponent, Sz, ruled by the Co 3d electronic states lying 
within the bulk band gap. This Sz component is ^100% 
polarized, whereas the polarization of the surface heli- 
cal spin-texture of the massive Dirac cone is drastically 
reduced with respect to pristine Bi2Se3. The impurity 
introduces a net spin with MAE of 6 meV per adsorbed 
Co atom, favoring the magnetic moment to be aligned 
perpendicular to the TI surface. These results have a 
direct connection with information processing and spin- 
tronic applications using topological devices. 
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